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The pairwise velocity probability density function in 
models with local primordial non-Gaussianity 
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ABSTRACT 

We study how primordial non-Gaussianity affects the pairwise velocity probability 
density function (PDF) using an analytical model and cosmological iV-body simula- 
tions. We adopt the local type non-Gaussian models characterized by /„/, and examine 
both the linear velocity difference PDF and the linear pairwise velocity PDF. We show 
explicitly how /„; induces correlations between originally independent velocities along 
the parallel and the perpendicular to the line of separation directions. We compare 
the model results with measurements from A^-body simulations of the non-Gaussian 
models. Linear theory fails to predict the PDF in the /„/ models. Therefore we develop 
an analytic model based on the Zeldovich approximation to describe the evolution of 
the velocity PDF. Our analytical model and simulation results show remarkably good 
agreement in both the parallel and the perpendicular directions for the PDF profiles, 
as well as the change in the PDF due to primordial non-Gaussianity. The agreement is 
particularly good for relatively small separations (< 10/i _1 Mpc). The inclusion of the 
evolution of the velocity PDF is important to obtain a good description on the signa- 
ture of primordial non-Gaussianity in the PDF. Our model provides the foundation 
to constrain /„/ using the peculiar velocity in future surveys. 

Key words: methods: analytical - dark matter - large scale structure of the universe 



1 INTRODUCTION 



Primordial non-Gaussianity has attracted much attention owing to it s ability to distinguish inflationary models fe.g. jBuchbinder et al 



20USl : lKhourv fc Piazzahqoll Silvestri fc Troddenl2008l:IWandsl20ld. and references therei n). Sev eral cosmol ogical probes using 



the CMB (Komatsu et al. 



201d:lHikag e et aljl2008l : lYadav fc Wand clt 2008; McEwen et al.ll2008l:lRossi et al.ll20 09: Smid t et al 



2010 |: lKomatsull2010l:l Ro^i*^tal*02Olo]) and the large-scale structures in the Universe iK^vjmiaetidT^ 



200C JScoccimarro et al.ll2004l;ISefusatti fe Komatsul2007l;llzumi fc Sodal2007l;lLo Verde et al.l2008l;lDalal et al.l2008l;lMatarrese fc Verdel 



2008 



Carbone et al 



Grossi et al 



20091 : 
20ld: 



2008 



20081: lAfshordi fc, Tollevl |200S|; Isiosar et all [20081: iMcDonaldl 20081: iTaruva et al 



Kamionkowski et al 



Lam et al. 2009, 201 



10: 



Tseliakhovich et al] 201 



20091: IPesiacaues et al. 
~2009l:ICunha et al.l 



km. 

2010; 1 



Smith et alj|2010l ) have been considered 



Nishimichi et al 
d lXia et al"ff2010l : IChongchitnan fc Silkl2010l 



Pillcpi ch et al 
Sartori^ e^alJI^Ol 



2008 



Lam fc Sheth 



2008 



2009 



Tsiosarl 120091: 
jGrossi_et_alJ_ 



DcsiacQucs 



Scliak 2010a; Liguori et al 
Giannantonio fc Porcianill20ld : l]VLadan^taITl201ol : 



Here, we introduce a new large-scale structure probe of primordial non-Gaussianity - the pairwise velocity PDF. Current 
large-scale structure probes of primordial non-Gaussianity focus on the clustering of halos, halo/void abundances, bispectrum, 
and the PDF of dark matter field - all of them are related to the change of the density field due to primordial non-Gaussianity. 
It is important to note that both the initial density field and the initial velocity field are generated from the primordial 
perturbation, and that their linear relation is described by the continuity equation. Hence, one may naively expect that 
primordial non-Gaussianity affect the velocity field as well as the density field. While various measurements associated with 
the change of densit y field have been extensively studied, there are few studies on the effect of primordial non-Gaussianity 
on the velocity field. Scherrer 1 19921 ); Catelan fc Scherrer (1995) discussed the effect of primordial non-Gaussianity on the 
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distributions of linear velocity fields. In par ticular Catelan fc Scherrer ( 19951 ) studied the distribution of the parallel component 
of the pairwise velocity using linear theory. Lam et al.l ( 201fJ ) discussed the effect of primordial non-Gaussianity on the redshift 
space distortion. They made use of the ellipsoidal collapse model to derive the modification in the Kaiser factor relating the 
real a nd redshi f t spac e power spectra, without discussing the modification in the velocity field. During the preparation of this 
work, ISchmidtl l|201ol ) studied the primordial non-Gaussianity signature in the peculiar velocities of density peaks using linear 
theory. Our present study has two important improvements over these previous works: first we show that primordial non- 
Gaussianity induces a correlation betweeen velocities in the parallel and the perpendicular to the line of separation directions. 
We show explicitly how this correlation, which is absent in ISchmidt (2010), modifies significantly the linear velocity PDF. 
Secondly, we show that the linear theory does not provide a good description of the signature of /„; in the pairwise velocity 
PDF even at separation as big as 50 h~ Mpc. We use cosmological TV-body simulations to show this. We thus develop an 
analytic model to describe the evolution of the velocity PDF. Our model is based on the Zeldovich Approximation. We 

illustrate the improvement in the iV-body measurement comparisons. 

The peculiar velocity field has been invest i gated as a probe of cosmology for both dark matter jGorski 1988; Sct o fc Yokovama 



19981; ISheth et al.ll200ll;lKuwabara et ajj|2002l ; Iscoccimaxr5l20o4) and biased tracers dsheth fc Diaferiolbooi 



Sheth et al 



Hamana et al.ll2003l : ISheth fc Zehavill2009l ). Other studies (|Yoshida et al.1l200ll : |Peelll2006l : iBhattacharva fc Kosowskv 



2001 



2007 



200 



3) discuss the possibility of constraining dark energy via the kinetic Sunyaev-Zeldovich effect. In this study we focus on 



the peculiar velocity field of the dark matter field. The effect of primordial non-Gaussianity on the velocity of biased tracers 
will be discussed in future work. 

We first describe how the linear velocity PDF changes due to primordial non-Gaussianity in section [2] Throughout the 
present paper, we work with models with non-vanishing primordial bispectrum; we will use the local f n i type primordial 
non-Gaussianity to illustrate the calculations. The Bardeen potential $ in the local /„; model is 

<6 = </> + / ni (0 2 -<0 2 », (1) 

where <j> is a Gaussian potential field and /„; is the nonlinear quadratic parameter. The above non-Gaussian correction is 
defined at z = zcmb for this study. It has been suggested that mass weighting is important for the peculiar velocities 
1 Scoccimarrol 2004; Sheth fc Zehavi 20091 ) . We will discuss the effect of primordial non-Gaussianity on both the uniform 
weighted as well as the mass weighted linear velocity PDF. In section [3] we describe the analytical model to approximate the 
evolution of the pairwise velocity PDF. The theoretical predictions of both the linear theory and the analytical model are 
compared with measurements from iV-body simulations. We conclude our findings in section [3] 



2 LINEAR VELOCITY PDF 
2.1 Preliminary 

The linear overdensity S(k) and the Bardeen potential <E>(fe) is related by the Poisson equation 

S(k, z) = D(z)k 2 M(k)$(k), (2) 

where M{k) — 2c 2 T(k)/3Q m H 2 and T(k) is the matter transfer function (note we do not include the k 2 factor in M(k)). The 
continuity relates the linear overdensity and the peculiar velocity u(k, z): 

j(fe) + 0(fc) = O, (3) 

where 6(x) = V ■ u(x) is the divergence of the velocity field. The velo city field is described so lely by its divegerence since its 
vorticity decays due to the expansion of the universe (for example see Bernardeau et al. 20021 ). Hence 



Uj(k) = W{z)k j M(k)^{k), (4) 

where i 2 — — 1 and the subscript j denotes the coordinate of the peculiar velocity. In this study we will be interested in the 
relative velocities in the parallel (||) and perpendicular (_L a and _Lf,) to the line of separation for two particles separated by 
some distance r. 

Connected moments higher than the second order vanish when the Bardeen potential is Gaussian. When the primordial 
perturbation is non-Gaussian, the leading order of the non-vanishing connected moment (higher than the second order) 
depends on the particular model of primordial non-Gaussianity. Most studies in the literature concern with models with a 
non- vanishing bispectrum (including the local f„i and the equilater al triangle type f n i)- Some ot her studies investigate models 



with a non- vanishing trispectrum (the g„i model, see for example iDesiacques fc Selialdl2010bl ). While this work focuses on 
primordial non-Gaussianity with a leadingly non-vanishing bispectrum, it can be extended to study models where the leading 
order of the non-vanishing connected moment is higher than the third order. In what follows we will use the local f n i model 
to demonstrate how the non-vanishing third order connected moment modifies the pairwise velocity PDF. The calculation 
can also apply to models with other types of primordial bispectrum. 

We denote the peculiar velocity at position x by u(x) and the relative velocity by v(r) = u(x) — u(x') = u — u' where 
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x — x' — r. We also denote the velocity difference PDF p(v; r) as the PDF of the peculiar velocity difference at two random 
positions separated by r; the pairwise velocity PDF q(v; r) as the PDF of the difference of the peculiar velocities of two tested 
particles separate d by r. The latter PDF is the pair weighted version of the former one and they are in general not equivalent 
(see, for example Scoccimarrol 2004). For simplicity we sometimes call the former PDF the uniform weighted PDF and the 
latter one the mass/pair weighted PDF. 

The numerical integra tion of the third order connected moments are carried out by Monte-Carlo integration using the 
numerical package CUBA (Hahn 2005). 



2.2 Linear Velocity Difference PDF 

2.2.1 Case: f nt =0 

The linear velocity difference PDF of two random positions separated by a distance r when /„; = is given by the multivariate 
normal distribution: 



Po(v;r) 



1 



exp 



1 T a-1 

-—v A v 



(5) 



(27T)S/VM 

where v = (vu,v± a ,v±.) in which Vii corresponds to the relative velocity parallel to the line of separation and v_i_ a and v± b 
are the two velocity differences perpendicular to the line of separation. Here A is the covariance matrix. The above expression 
simplifies since there is no correlation between different components of v. Hence the r.h.s becomes a product of three univariate 
normal distributions: 



Po(«n , v± a , vx b ; r) = p («|| ; r)p (vj_ a ; r)p (vj_ b ; r), 
where the variances of the univariate normal distributions are and (t>j_) = {vj_ ) = (vj_ ) respectively, and 



(tfi) = ^VoJ dfcP*(fc)fc 4 M 2 (fc) 
(° 2 ^ = j dkP*(k)k 4 M 2 (k) 



3jo(fcr)+6 



j i(kr) 
kr 



kv 



(6) 

(7) 
(8) 



Here jo and ji is the spherical bessel function, Do = dD /At is the time derivative of the linear growth factor D at z — and 
P<s>(k) ~ P<t>{k) is the Bardeen potential power spectrum. Iostropy means one can transform the rectangular coordinate in the 
plane perpendicular to the line of separation into the polar coordinate and write equation ((6| as 

Po(v\\ , v± a , v± b ; r)dv^dv± a dv± b = 2irv±po(v\\ ; r)po(v± ; r)dv^dv± , (9) 



where v\_ — v\ a + v ± b - Notice that our notations (v± ■ ■ ■)(= (v± a ■ ■ ■) = (v± b ■ ■ ■)) and v±(— \J v \ a + v ± b ) 



are not equivalent. 



2.2.2 Case: /„; / 

When f n i 7^ 0, the primordial bispectrum is non-zero and its functional form in the local type is 

B (ki, fca, few) = 2f nl [P(k 1 )P(k 2 ) + cyclic] + 0{fli). (10) 

Connected moments higher than the second order are non-vanishing and contribute to the velocity difference PDF. The 
leadingly non- vanishing connected moments are 
(v'n) ior>3 



12A1 



d :i fci 



2f nt (2tt) 6 

(y^l) _ (v\\vi a ) _ ( v \\ v l b ) 

2/n! 



cos ^2 5:0 
2 



d J fe 2 P(fc 1 )P(fc 2 )M(fc 1 )M(fc 2 )M(fc 12 )fc 1 ||A; 2 ||fc 12 ||[sin(fe 12 ||r) -2sin(fe 2 ||r)] 



(11) 



d 3 fci 



d 3 fe 2 P(fci)P(fc 2 )M(fci)M(fc 2 )M(fci 2 )[fc 1± fc2±fci 2 || sin(fci 2 ||r) 



2/ni 2/ ni (2tt) 6 j Jcos M 2>0 

- 2fc 1± fc 2 ||fci 2 ± sin(fc 2 ||r) + 2k 1± k 2 \\k 12 ± sin(fc 12 ||r) - 2fci|| fc 2± fci 2 ± sin(fe 2 ||r) - 2k 1± k 2 ±k 12 \\ sin(fe 2 ||r)], (12) 

where k\ 2 — fci+fc 2 |, ki\\ = ki cos /ij, ki± = ki sin fii cos cj>i, and fe 12 i| =fei||+fe 2 || (similarly for ±). Note that connected moments 
involving odd power of vj_ vanish irrespective of the value of /„(. The non- vanishing component (uiit)j_) results in couplings 
between velocities in the parallel and the perpendicular to the line of separation directions. As a result the linear multivariate 
PDF can no longer be written as a product of three independent univariate PDFs. We use the tri-variate Edgeworth expansion 
to approximate the linear uniform weighted velocity difference PDF for /„( 7^ (see appendix for derivation): 

p(v\\,vj_ a ,v± b ;fni,r) = Po(v\\,vj_ a ,vj_ b ;r)[l + 0:300/1300 + "120(^120 + ^102)], (13) 

where 



a 300 



1 Hi) 

6 {vfy/ 2 ' 



1 

2 {vf,)V 2 {t>l}' 
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J i i i I i i i i i i 

2 4 6 

ln(r) 



Figure 1. Scale dependence of parameters: «3oo (thick solid black curve), ai20 (thick dashed red curve),— /3ioo (red sketeal trianglar 
symbols), /3200 (thick cyan short-long-dashed curve), —7100 (green starred symbols), 7300 (thick green dot-dashed curve), 7120 (thick 
blue dotted curve), —7400 (blue pentagons), —7220 (cyan hexagons), 7500 (thin magenta short-long-dashed curve), and 7320 thin (orange 
dot-dashed curve). /„; = —100 for those parameters apply. 0300, cti20i — 7100; and — P100 are shifted vertically by -2 for clarity. 



and hijk = Hi(y\\)Hj(yx a )Hk(yx h ) is the product of Hermite polynomials of different orders. In particular, 



,.3 



^300=^-3^1 = ^5^-3^3^ (14) 



i(^- l ) = (JW^- 1 ) 



^ 2 = ^„K b -l) = ^^-lj. (16) 

Note that ISchmidtJ [j2010h marginalized velocities perpendicular to the line of separation, essentially setting cv.120 = 0. 

Figure [T] shows the scale dependence of 0300 (thick solid black curve) and Q120 (thick dashed red curve) for /„; = —100. 
While the numerical value is small compared to unity, the effect of non-zero /„; becomes significant for big \vn \ or \v±\. 
Furthermore CV300 and 0120 are in the same order of magnitude, hence both terms have to be included in the computation of 
the linear peculiar velocity difference PDF. 

Figure [2] shows the ratios of the linear peculiar velocity difference PDF for f n i = 100 to the corresponding Gaussian PDF 
at r = 8 /i^Mpc. The axes labeled vll indicates the velocity in the parallel to the line of separation direction and vL is the 



magnitude of the velocity perpendicular to the line of separation (recall vL = \J v \ a + v j_ b )- The upper panel shows both the 
contour and color maps over a wide range of velocity. The effect of primordial non-Gaussianity is most significant at extreme 
velocities. Both the contour and the color maps show variations in the ratio for different values of v± at fixed v\\, indicating 
that the modification due to f n i is degenerated in the plane of (vll, vL). While the upper panel shows that the ratio can be as 
big as 4, the first order Edgeworth expansion breaks down in such rare cases (about 10-ct level) - it is evident by the negative 
probability at the other end of the plot. 
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Figure 2. Ratio of linear velocity difference PDF p/po in equation (1 1 3 I t for r = 8 h~ 1 Mpc and / n ; = 100. vll indicates velocity in the 
parallel to the line of separation direction while vl_ is the magnitude of the velocity perpendicular to the line of separation. Upper panel 
shows the ratio for a wide range of velocity. The lower left panel zooms into regions for around 5-<r as indicated by the outer axes. The 
lower right panel shows the same regions as the lower left panel, but explicitly set the parameter C1120 = as previous study did. 

The lower panels of Figure [2] show the color maps of the same PDF ratio on a smaller velocity range to increase the 
dynamic range of the effect of /„; on the linear velocity difference PDF. The left panel shows the modification obtained in 
this study (the correction term in equation (|13[) 1 while the right panel shows the result for setting (wyDx) — 0. The outer axes 
of the lower panels show the tr-level in corresponding directions. At 5-a level (vll direction) the change is about 30%, and the 
bigger the <r-level in the perpendicular direction the greater the modification is. The linear theory suggests that for /„; > 
the infalling probability (negative velocity) is enhanced while the outgoing probability (positive velocity) is decreased. This 
trend applies to other scales we looked at (from 4 to 100 ft~ 1 Mpc). While it is generally believed that linear theory applies 
for large scales, we will show in the next section that linear theory fails to describe the change in the velocity PDF due to f n i 
even at separation as big as 50 ft _1 Mpc. 

2.3 Linear Pairwise Velocity PDF 

It is widely accepted that, while the linear theory prediction is a good approximation for large scales (> 20 /i _1 Mpc), the 
linear theory velocity correlations of massive halos (in the parallel to the line of separation direction) is not consistent with 
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iV-body measurements (|Croft fc Efstathioulfliiih . Isheth fc Zehavil (|2009h pointed out that the velocity correlations of biased 
tracers, in both parallel and perpendicular to the line of separation directions, can in fact be reasonably described by linear 
theory when proper pair-weighting is included. While our current study focuses on the difference of the peculiar velocities of 
unbiased tracer, it is useful to examine whether the pair weighting would improve the linear theory prediction. It is reasonable 
to include this mass weighting as the pairwise velocity PDF measured in iV-body simulations is obtained by counting the 
relative velocities of simulated particles. As a result regions with more particles (overdense regions) have a bigger weight than 
regions with less particles (underdense regions). 



2.3.1 Case:f„i=0 

Scoccimarrol (2004) discussed how to generalize the linear velocity difference PDF when f n i — to include the mass weighting. 



The resulting linear pairwise velocity PDF qo is related to the associated no weighting velocity difference PDF po by 



[l + £(r)l 



1 + £(r) + ftiooAoo + ^200^200, 



where 8 



1 po(u h u ±a ,u ±b ;r) 

: S(x), 8' = 8(x + r), £ = £(r) = {85') is the matter correlation function at separation r, and 



0100 = 



(vnS) 



(W|l<5) <V||<5'} 



2\l/2 ' 



0200 = 



{v\\6){v\\5') 



27T 2 ' 



Do / dkP ii (k)k 4 M 2 (k)[kj 1 (kr)}. 



(17) 



(18) 



(19) 



The r.h.s of equation ()17p does not depend v±, hence qo can still be written as a product of three univariate PDFs - the 
correction term on the r.h.s of equation (|17|) only modifies the univariate PDF of vu . The scale dependence of the parameters 
— 0ioo (red sketeal triangular symbols) and 0200 (thick cyan short-long-dashed curve) are shown in Figure [T] their magnitudes 
are big compared to others shown in the same figure since they are 'Gaussian parameters' that do not depend on f n i. 



2.3.2 Case:f n i^0 

When /„; 7^ additional terms contribute to the mass weighted linear velocity difference PDF. To the first order of f n i the 
expression is (see derivation in appendix) 

[1 + C( r )] g ^ V ±a ' Xb '^ n '' ^ = 1 + £(r) + ftiocmoo + /1200/3200 + ^3007300 + (/1120 + ^102)7120 

+ /l4007400 + (^220 + /l202)7220 + /l5007500 + (^320 + /l302)7320, (20) 

where 

7ioo = Pioo + pyiTS 7300 = asoo(l + t(r)) + — -^y 7120 = Qi2o(l + C(r)) + — -^y 

0100 a a 

7400 = Q300 ~ 7220 = «120P100 7500 = Q300P200 

7320 = O120/3200 h,A00 = V\\ ~ + 3 hb00 = V\\ — lOfc'l + 15f || , 

and for the local /„; type primordial non-Gaussianity 

^rl^ = -75% / d3fel / d 3 fe 2 P(fc 1 )P(fc 2 )Af(fc 1 )Af(fc 2 )M(fc 12 )[fc 1 ||fcifc 12 sin(*i 2 nr) + fc 2 fc 2 2 fci 2 || sin(fc 2 ||r) 

2,/ni (27T) J Jcos M2 >0 

- fci||fc 2 fc? 2 sin(fc 2 ||r)] (21) 
■hrr 1 = 77T^ / d3fc i / d 3 fe 2 P(fci)P(fc 2 )M(fci)M(fe 2 )M(fei 2 )[2fc 1 ||A; 2 2 fc 12 || cos(fe 2 ||r) - fci||fc 2 ||fc 2 2 cos(fc 12 ||r) 

2/ni (27T)° J Jcos M2 >0 

+ 2fc 1 ||fc 2 ||fc 12 cos(fc 2 ||r) - 2fc 1 ||fc 2 fci 2 || cos(&i||r) - 2fe 1 ||fc 2 ! fei2|| cos(k 12 \\r) + 2fc 1 ||fc 2 fc 12 || - fci|| fc 2 || fc 2 2 ] (22) 
= TT^E [ d3fel / d 3 k 2 P(k 1 )P{k2)M{k 1 )M{k 2 )M(k 1 2)[2k 1± k 2 2 k 1 2± cos{k n r) - k 1± k 2± k 2 12 cos(k m r) 

+ 2k 1 _ L k 2 ±k 2 2 cos(fc 2 ||r) - 2fcixfc|fei 2 ± cos(fci||r) - 2k 1± k 2 ki 2 ± cos(fc 12 ||r) + 2fci||fcf fc 12 || - fci||fc 2 ||fc 2 2 ]. (23) 

The scale dependences of jijk are shown in figureQ] As in the case of uniform weighted linear PDF, the linear pairwise velocity 
PDF can no longer be written as a product of three independent PDFs when /„; ^ due to the non- vanishing terms (v|| t)j_ ) 
and (t)j_<5'). In addition, when one sets 8 = 8' = 0, the above expression recovers the uniform weighted linear PDF for f n i 7^ 
(equation 1)130). 

Figure [3] shows the modification factor of the linear pairwise velocity PDF due to fni (q/qo in equations (|17|) and ([20])). As 
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Figure 3. Similar plot to figure[2]but for the ratio of the linear pairwise velocity PDF (q/qo in equations II17I I and 1 1201) ), The lower right 
panel sets a\20 = 7220 = 7320 = 0>x^') = °- 



in Figure [2] the upper panel shows the ratio in a wide range of velocity, including regions with negative probability indicating 
the breakdown of the first order Edgeworth expansion approximation. The lower left panel enlarges the color map on regions 
about 5-(j level and the outer axes label the a-level. The lower right panel shows the ratio if one explicitly sets all third order 
moments involving v±_ to zero (i.e. «i2o = 7220 = 7320 = <oi<5') = 0). The effect of f ni in the mass weighted linear pairwise 
PDF has the same trend as the uniform weighted PDF: the effect of f n i is degenerated in the parallel and the perpendicular 
to the line of separation directions; and the general trend observed in uniform weighted linear PDF still applies - for /„; > 
the probability of finding infalling pairs increases while the probability of finding pairs moving apart decreases. The inclusion 
of the mass weighted quantites 7^ strengthens the effect of f n i and it is most significant in the decrement in probability of 
outgoing pairs. 



3 EVOLVED VELOCITY PDF 

Given the linear velocity PDF we now describe a model for the evolution of the PDF. We adopt the Zeldovich Approximation 
which assumes the comoving velocity remains unchanged. At some initial redshift Zi, if two particles separated by n have a 
relative velocity (vl , Wj_ , v l ± ). Then at a later redshift z = zq (<C Zi) the distances traveled are Da/Di(vpV t La ,v t ±b ) (note 
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that it assumes Do S> -Dj). Hence the separation of the two particles becomes 

2 / j-, \ 2 



2 

r 



r! + ^"J Ha 1 +u±b } ' (24) 



The evolved relative velocities (with respect to the updated position) are: 

D ( nvl Do ,2 



I 1 2 2,2 / A) A 2 fr> a \ 

Kl = V ±a + V ± b = I ) ( 26 ) 

'2 -2 2 -2 

where t> 1 = vl + iq^ + ^l,, ■ The evolved velocity difference PDF is therefore 

p(Vj|,Vl;fl) = j drid^d«i s di;i b ^p(t>{S,i;i s ,i;i b ;r)fo(r-/i)fo(«||-Vi|)fo(«x-Vl), (27) 

(see, for example. ISeto fc Yokovamalll998h where the dirac-delta functions use equations (|24[) . (|25|) . and (|26|) to map the initial 
quantities to the evolved one. In the following we apply the above formula to compute the evolved velocity difference PDF 
from equations ([6]) and (|13[) and compare the results to measurements from TV-body simulations. We set Zi and zo to be 1100 
and 0.5 respectively, justifying the approximation Do S> Di. 



3.1 Comparisons to TV-body measurements 

We measured the pairwise velocity PDF from a set of TV-body simulations described in Nishimichi et al. ( 20091 ). The simulation 



adopts the WMAP 5-years ACDM best fit parameters (fi m , fi A , M b , h, cr 8 , n„) = (0.279,0.721,0.046,0.701,0.817,0.96). The 
simulation was performed in a box of 2000 7i _1 Mpc on a side, containing 512 3 particles whose mass is 4.6 x 10 12 /i _1 M . The 
measurements were made from the simulation output at z = 0.5. To account for the discrete nature of TV-body measurements, 
separation r in the following refers to two particles having a separation between (r — 2,r + 2) in the unit of h~ 1 Mpc. For 
the theoretical model we apply the same separation selection in the dirac-delta function of r. The comparisons of the velocity 
PDF of vu and v± are presented in the next two subsections. We will then discuss the implications of the comparisons in the 
next section. 



3.1.1 Parallel to the line of separation: p(vu;r, f n i) 

Figure|4]shows the PDF comparisons of the various analytical predictions to the iV-body measurements for the relative velocity 
parallel to the line of separation for 4 different separations (4, 8, 12, 50 /i _1 Mpc). The upper panel in each subfigure shows 
the PDF profile for f n i = 0: solid symbols are the TV-body measurements, cyan dotted curves are the analytical predictions of 
equation (|27|) (marginalized over v±), green dot-long dashed (mass weighted) and magenta long dashed (uniform weighted) 
curves are the linear theory predictions respectively. 

The predictions from our analytical evolution model give good matches to the TV-body measurements at all scales we 
investigated, except at large outgoing velocities (i.e. Vu ^> 0). The agreements with the TV-body measurements are better in 
small separations (4 and 8 /i _1 Mpc) than large separation (12 /i -1 Mpc). The analytical predictions of our model are able to 
describe the high infalling velocity regimes (vu < —1000). This is very encouraging since one may expect that high velocity 
regions are not described by our simple model. The analytical model also predicts knee-shape transitions at around 300, 600, 
and 900km/s for r = 4, 8, 12 ft _1 Mpc respectively. Similar but less significant knee-shape changes can be seen from the 
TV-body measurements around the corresponding velocities. In the high outgoing velocity regions the match is not as good as 
the infalling regimes. We believe the disagreement in the outgoing velocity regime as compared to the infalling regime is due 
to non-linear evolution of velocity that is not described by our current model - originally infalling pairs eventually become 
outgoing (see equation (J25J ) and non- linear interactions between particles in close contact are not described by the Zeldovich 
Approximation. 

On the other hand, the linear theory predictions do not fare as well as the analytical model's. At small separations the 
predictions of the linear velocity difference PDF do not match the TV-body measurements - the magneta curves miss both the 
extrema and the peaks of the PDF at separation r = 4, 8, and 12 h~ Mpc. This is consistent with previous studies that the 
linear theory uniform weighted velocity PDF does not agree with the velocity correlation in the parallel to the line of separation 
direction. On the contrary the mass weighted linear predictions provide re asonable matches to th e TV-body measurements 
near the peaks of the PDF at different separations. It confirms the finding of Sheth fc Zeha vi (2009): the velocity correlatL 



tion 



can be described by the linear theory when the mass weighting is taken into account. The matching of the peaks of the PDF 
guarantees the velocity correlations, which is equivalent to the expected value of the pairwise velocity PDF, would roughly 
agree. However the mass weighted linear predictions does not match the TV-body measurements when | is of the order of 
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Figure 4. Pairwise velocity (parallel to line of separation) PDF at different separations. Negative velocity corresponds to particles 
infalling while positive velocity indicates particles moving away from each other. Upper panels show the pairwise velocity PDF when 
f n l = 0: solid symbols are measurements from the TV-body simulation, cyan dotted curves are the velocity difference PDF from the 
evolution model based on the Zcldovich approximation (equation JUjl), green dot-long dashed curves are the linear pairwise velocity (go 
in equation 117> ). and magenta long dashed curves are the linear velocity difference PDF (equation (|6]l). Lower panels show the ratios 
of the PDF for /„; ± 100 to the associated PDF for /„; = 0: crosses (/„; = 100) and solid squares (/„; = —100), cyan dotted curves 
{fnl = 100) and blue dot-short dashed curves (/„; = —100) are predictions of equation H27I I. green dot-long dashed curves are the ratios 
for the linear pairwise PDF (f n i = 100), and magenta long dashed curves are the ratios for the linear velocity difference PDF (/„; = 100). 
Red short- long-dashed curves show the prediction from equation (1271 . neglecting third order moments involving velocity perpendicular 
to line of separation (see the lower right panel of figure . 



a few hundreds km/s. The disagreement is getting worse as the separation decreases. At large separation (50 /i _1 Mpc), the 
two linear theory predictions are very similar and they describe the JV-body measurements well. They are very closed to the 
prediction of the analytical model and the match to the measurements is only slightly worse at high velocity regimes. 

The lower panels of figure [4] show the ratios (substracted by unity) of the pairwise velocity PDFs for /„; — ±100 to the 
associated PDFs for f n i — at different separations. Measurements from the iV-body simulations are represented by crosses 
(fni = 100) and squares (/„; = —100). Other curves in the lower panels show the PDF ratios of different analytical predictions: 
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cyan dotted (/ n ; = 100) and blue dot-short dashed (/ n ; = —100) curves are the predictions of the evolution model based on 
the Zeldovich Approximation (equation l|27fl); red short-long-dashed curves show the similar predictions for / n j = 100 but 
setting («|it)j_) = (see the lower right panel of figure [2] for the corresponding linear theory comparison at r = 8 ft _1 Mpc). 
The other two curves are the linear theory predicted ratios for f n i — 100: the magneta long dashed curves are the uniform 
weighted velocity predictions and the green dot-long dashed curves are the mass weighted velocity predicted ratios. 

For the range of separations we investigated the change in the pairwise velocity PDF due to /„; = ±100 is at most 5% 
in the infalling velocity regime and 10% in the outgoing velocity regime. In contrary to the linear theory predictions (see 
figures [2] and [3] and the magenta long dashed as well as the green dot-long dashed curves), the effect of positive f n i enhances 
the probability of having pairs in both the infalling and the outgoing high velocity ends. This is true for all the separations 
we look at, from 4 to 50 /i~ 1 Mpc - scales in which the linear theory is believed to be valid. We also check using a simulation 
with a smaller box (1 h~ Gpc on a side) to estimate the box size effect: shrinking the box volume to one-eighth of its size 
changes the PDF ratios by at most 2% at vu = 2000 km/s. 

The predictions of our evolution model, regardless of whether one explicitly set (vn t)^) = , match the iV-body measure- 
ments reasonably well when r = 4 /i _1 Mpc. The predicted ratios from the analytic model match the measurements within 
2% across the velocity range we look at. It is remarkable to have such agreement at velocity as high as ~ 2000km/s. The 
predicted ratios also match the transition of the PDF from the infalling to the outgoing regime: for /„; = 100 it gradually 
decreases from positive to negative, then changes direction and become positive again. The same transition is also observed 
at larger separations and the model's predictions are able to match the ./V-body measurements. In contrast the ratios of the 
linear theory predictions (both the uniform and the mass weighted) fail to match the iV-body measurements in both infalling 
and outgoing velocity regimes. While the linear theory predicted ratios still have the same trend as the iV-body measurements 
in the infalling velocity regime (but the differences between the predictions and the measurements are large), the predicted 
ratios from linear theory are opposite to the measurements in the outgoing velocity regime (the same is also observed at larger 
separations). Hence the linear theory cannot be used to predict the effect of /„; on the pairwise velocity PDF. 

At larger separations the agreements in the ratios of the PDF between our analytical model predictions and the N- 
body measurements are less impressive: while the predicted ratios still agree well with the measurements in the infalling 
velocity regime, the differences in the outgoing velocity regime are getting bigger as the separation increases. Nonetheless 
the predictions still qualitatively match the measurements at r = 8 and 12 /t _1 Mpc. At r — 50 ft _1 Mpc, on the other hand, 
the analytical predictions converge to the linear theory predictions and show opposite trends to the measurements. In the 
outgoing velocity region, there are some interesting transitions in the analytical model predictions: at r = 12 ft _1 Mpc and 
fni = 100, the predicted ratios (the cyan as well as red curves) first show a decrement at small v« . They then turn around 
near v\\ = 800km/s and keep increasing for larger Vu. Such transitions are not seen in the TV-body measurements. 

We will leave the discussion on the implications for the discrepancies between the linear theory predictions and the 
iV-body measurements at all separations as well as the mismatch in the PDF ratios between our analytical models and the 
measurements at large separations in the next section. 

3.1.2 Perpendicular to the line of separation: p(v±;r, f n i) 

Figure [5] shows the comparisons of the PDFs of the pairwise velocity perpendicular to the line of separation for different 
separations. The upper panels show the PDF at z = 0.5 for /„; = and the lower panels show the ratios of the PDF for 
fni — ±100 to the associated PDF for f n \ — 0. The symbol and curve labelings are the same as in figure [4] but we do not 
show the PDF ratio predictions of the linear theory in the lower panels since the linear theory predicts no effect of f n i in the 
marginalized PDF, that is p(v± ; r, f„i)/p(v± ; r, /„; = 0) = 1. 

The predictions of the PDF profile from our analytical model (equation (|27p ) match the iV-body measurements reasonably 
well, except around v± = at r = 50 /i _1 Mpc. The predictions agree with the measurements in a wide range of velocity 
(200 — 1500km/s), as well as the dips near v± = 0. At higher ui the model predictions always underestimate the PDF. The 
linear theory predictions, on the other hand, fail to match the measured PDF for r = 4, 8, and 12 /i _1 Mpc. At r = 50 fo _1 Mpc 
the linear theory prediction agrees with the iV-body measurements for ui = 300 — 1200km/s, but it fails to predict the dip 
near v± — 0. Note that there is no difference between the mass weighted and the uniform weighted linear theory predictions 
since the correction term due to mass weighting does not depend on v± (equation (|17[) ). 

The effect of primordial non-Gaussianity in the pairwise velocity PDF in the perpendicular to the line of separation 
direction is not as strong as the signature in the parallel to the line of separation direction: the change in the PDF is top at 
7% at v± = 2000km/s. The ratio first shows a decrement (increment) at small v± for /„; = 100 ( — 100). It then gradually 
switches direction and at high v± show an increment (decrement). This gradual switch from decrement to increment is 
most significant at small separation. Our analytical model is able to predict the change in the PDF due to primordial non- 
Gaussianity qualitatively. Its prediction matches the gradual change of the PDF ratios (from decrement to enhancement for 
fni = 100, opposite for f n i = —100) and provides a rough approximation on the velocity where the crossing across unity 
occurs. However the predictions always underestimate the PDF ratios in low v± regime. At higher v± the model prediction on 
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Figure 5. Pairwise velocity (perpendicular to line of separation) PDF at different separations. Upper panels show the PDF profile when 
fnl = : solid symbols are measurements from the ./V-body simulation, cyan dotted curves are the prediction from the evolution model 
based on the Zeldovich Approximation (equation 1270 ). green dot- long dashed curves are the linear pairwise velocity (go in equation JT7J), 
and magenta long dashed curves are the linear velocity difference PDF (equation {6}). Lower panels show the ratios of the PDF when 
fnl ± 100 to the associated PDF for /„; = 0: crosses (/„; = 100) and solid squares (/ n ; = —100) are measurements from TV-body 
simulation, cyan dotted curves (/„; = 100) and blue dot-short dashed curves (/„; = —100) arc predictions of equation I I27II . Red short- 
long-dashed curves show the prediction from equation {27} . neglecting third order moments involving velocity perpendicular to line of 
separation (see the lower right panel of figure [2} 



the PDF ratio matches the JV-body measurements very well at r = 8 as well as 12 /i -1 Mpc; the predicted ratio overestimates 
the change at r — 4 /i -1 Mpc but underestimates at r = 50 /i _1 Mpc. The analytical model predicts the effect of f n i becomes 
weaker at some fixed v± for larger r, the measurements from the iV-body simulations, on the contrary, remain roughly the 
same at various separations. The measured PDF ratios at high v± only depends on the scale of separation very weakly. In 
addition there are spikes near v± — in the theoretical predicted ratios for r = 12 and 50 /i _1 Mpc, but these are not seen in 
the iV-body simulations. 

The linear theory predictions fail to match the TV-body measurements at all separations - it predicts no change in the 
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pairwise velocity PDF in the perpendicular to the line of separation direction. The TV-body measurements, on the other hand, 
has several percent-level change in the PDF throughout the velocity range we look at. 



4 DISCUSSIONS 

We have studied how primoridal non-Gaussianity affects the probability density function of the pairwise velocity. Our study 
complements earlier studies on the signatures of primordial non-Gaussianity on other large-scale structure probes which utilize 
the modification in the density field (for example the PDF of density field, halo/void abundances, scale dependent halo bias). 
We have shown that primordial non-Gaussianity models with a non-vanishing bispectrum result in two non-vanishing third- 
order connected moments of the peculiar velocity and (ti||t>j_)). We adopt the local f n i type primordial non-Gaussianity 
as an example to compute the modification in the linear velocity difference PDF. Following suggestions by earlier studies, 
we investigate the effect of /„; on the linear pairwise velocity PDF. Primordial non-Gaussianity induces correlations between 
velocities in the parallel and the perpendicular to the line of separation directions. Hence both the uniform weighted and the 
mass weighted linear velocity PDF can no longer be written as a product of the corresponding univariate PDFs. Futhermore 
the change in the linear velocity PDF is degenerated in the plane. 

Next we have developed an analytical model to describe the evolution of the velocity PDF. It employs the Zeldovich 
Approximation which assumes the comoving velocity remains unchanged. We then compared the pairwise velocity PDF 
measured from a TV-body simulation with the predictions from our analytical model at various separations, ranging from 4 
to 50 ft _1 Mpc. Our model predictions agree well with the PDF profiles measured from the TV-body simulation at various 
separations when /„; = in both the parallel and in the perpendicular to the line of separation directions. 

Linear theory fails in matching the TV-body measurements: both the linear velocity difference PDF as well as the linear 
pairwise velocity PDF fail to predict the PDF profiles and the change in the PDF due to a non-zero f n i- The linear theory 
prediction does not take into account the movement of particles: pairs separated by r having a relative velocity (vu,v± a ,v±,) 
must have had different separations at earlier redshift. The relative velocities at earlier redshifts also change in general due 
to the difference in the particles' positions. In other words linear theoy is a good approximation only when the displacement 
of the particles is small compared to the separation between them. For example the displacement associated with a relative 
velocity v — 1000 km/s is approximately 10 fe _1 Mpc at z — 0.5 if the velocity remains unchanged. Therefore, only for large 
enough separations (r ^> 10/i~ 1 Mpc) the linear theory can approximate the PDF. The initial separation and the initial 
relative velocity determine the evolved separation and the pairwise velocity - equations (|24|l to (126 p give an example of such 
an evolution model using the Zeldovich approximation. In general the mapping from the initial parameters {r l , vl , v l ± } to the 
evolved parameters {r, v\\ ,v±} is non-linear and more than one set of {r l , , v l ± } would evolve to {r, wy , v± }. This analytical 
model provides results that match the measurements much better than the linear theory. 

Including evolution is also important in explaining the change in the PDF due to /„; . The ratios of the pairwise velocity 
PDF in the parallel to the line of separation direction can be used to reveal the distribution of pairs with different {r l , vl , v± } 
having evolved to the same {r,v«,Vj_}. Recall that the linear theory predicts enhancement in the probability occurs in high 
velocity infalling pairs and low velocity outgoing pairs when /„; > 0. On the other hand, the TV-body measurements always 
have increment in the probability for outgoing pairs. The analytical model provides an algorithm to explain this discrepancy: 
pairs initially infalling becomes outgoing when they cross each other on the projection on the line of separation. As a result 
outgoing pairs would be either originally infalling or outgoing pairs. The converse is not true because infalling pairs can only 
come from pairs that approach each other at earlier redshift. The originally low velocity outgoing pairs, having a higher 
probability for / n ; = 100, by itself cannot explain the observation from TV-body measurements as the change in the velocity 
difference PDF associated with these pairs is relative small (recall that the change in the PDF goes as w| — 3ii|| and v\_ — 1) 
and cannot account for the level of change in the PDF measured in the simulations. As a result the PDF ratio in the outgoing 
velocity regime has to be dominated by original infalling pairs. 

Our analytical model provides an approximation for the evolution of the velocity difference PDF and it successfully 
predicts the modification in the velocity PDF due to /„; qualitatively. The Zeldovich Approximation, on which the analytical 
model based, does not take into account several important features - it assumes the comoving velocity remains unchanged. 
This is a good approximation when the separation of the pair in consideration is big. This assumption may not hold when the 
pair crosses on the projection along the line of separation. Non-linear model is required to describe the evolution of the relative 
velocity in this cases. As a result the velocity PDF in the regime where cross-over dominates may not be well described by 
our analytical model. 

Another area our analytical model does not describe is the virial motion of particles residing inside dark matter halos. 
While it is not expected to have pairs residing within the same dark matter halo for the separation we considered in this 
study, virial motions would certainly affect the distribution of the pairwise velocity. 

Notice that our analytical model only describes the evolution of the velocity difference PDF - it is the evolved uniform 
weighted velocity PDF. In practice it is not what we measured from the TV-body simulation - it is a (evolved) pair weighted 
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velocity PDF. The comparisons with iV-body measurements show that the evolved velocity difference PDF matches the 
measurements better than the linear pairwise velocity PDF. Modeling of the non-linear pairwise velocity PDF requires the 
correlation between the non-linear density contrast and the evolved velocity field. These are not described in the our model 
and will be left for future study. 

It is quite unexpected that the agreement of the analytical model with the TV-body measurement is better at small 
separation than at large separation. One speculation for that is the simulation measurement is affected by the finite box 
size effect for big r - we did check the effect of finite box size by measuring the pairwise velocity PDF from a simulation 
ran in a box with 1 /i _1 Gpc on a side. The resulting PDF at r — 50 /i _1 Mpc agrees qualitatively with the corresponding 
PDF in the main text. Hence we believe the disagreement in the PDF ratio at r — 50 h~ Mpc is not due to the effect of 
finite box size. Another possibility is the discrepancy is due to the non-linear interactions - here we give two approaches to 
explain qualitatively how the non-linear interactions affect the velocity PDF. Both are related to the formation of non-linear 
structures and is not described by the Zeldovich Approximation. 

The Zeldovich Approximation assumes the comoving velocity is constant throughout the evolution and it is not a good 
approximation when non-linear structures (pancake, filament, and halos) form. If the Zeldovich Approximation were the 
correct model that completely described the evolution, the positive velocity in the evolved velocity difference PDF would 
be dominated by originally infalling pairs only when vu > Vo, where vo corresponds to the velocity that has a displacement 
approximately equal to the separation r. This scale also associates with the knee-shape transition in the PDF profiles in the 
upper panels of figure 3] For Vu < Vq the Zeldovich Approximation predicts originally outgoing pairs should dominate the 
PDF. Measurements from TV-body simulation, on the other hand, show no such transition. The knee-shape transition in the 
PDF profiles (upper panels) is also much weaker compared to the analytical predictions. Non-linear evolution models, such as 
the adhesion model, may be able to explain the difference between the measured and the predicted PDF - while the Zeldovich 
approximation predicts the velocity PDF is dominated by originally infalling pairs only when vu > Vo, the frictional dragging 
term in the adhesion model would damp the velocity, making wo < v'« —¥ V\\ < i>o. As a result the velocity PDF would be 
dominated by originally infalling pairs at a smaller vv, compared to the Zeldovich Approximation prediction. Notice that this 
effect will be most significant for big r as the dragging term is proportional to v. 

The halo model provides another explanation for this discrepancy: in the halo model all particles are residing in halos, 
each of them acquires its pecul iar velocity from two contributions - the virial motion inside the p arent halo as well as the bulk 
movement of the parent halo (jCoorav fe Shetblbooi ISheth fc DiaferiolboOll : Isheth et alJboOlh . Hence the relative velocity 
between pairs residing in different halos (this corresponds to the separations we studied in the current study) can be separated 
into two parts: the randomly-oriented velocity difference due to the corresponding virial motion and the relative velocity of 
the two parent halos. The dispersion of the virial motion depends on the mass of the parent halo: more massive halos have a 
larger dispersion; the bulk velocity of the pa rent halo depends only weakly on mass. The bulk velocity can be approximated 
by the linear theory ijCoorav fc Shethl l2002h . At relatively small separations, the two contributions are comparable; at big 
separations, on the other hand, the two-point correlation function is weak and the velocity is dominated by the virial motio n 
term. It is a well-known fact that there are more massive halos for models with f n i > (see, for example, Lam et al. 20091 ) . 



hence the velocity dispersion is bigger when / n j > 0. This effect is most significant in the pairwise velocity PDF for big r since 
at such separation the dominant velocity difference is the random-oriented virial velocity difference. The net result is having 
a higher probability when is big but lower probability when is small. Our current analytic model does not include 
the virial motion and as a result it is not able to describe the change in the velocity PDF for large r. Note that this picture 
naturally predicts the discrepancy should first show in the outgoing velocity regime - clustered halos attract each other and 
hence the relative bulk velocity is always infalling when projected on the line of separation of the pairs in consideration. 
Figure [6] shows the change in the PDF of peculiar velocity for models with f n i 0. The increase in velocity dispersion and the 
change in PDF at high v when /„; = 100 supports this halo model picture in explaining the discrepancy at large separation 
due to the increase in the magnitude of the virial motion in f n i{> 0) models. 

It is out of the scope of the current study to include these non-linear effect - this will be left for future study. Our current 
analytic model is able to predict the modification in the pairwise velocity PDF within a few percent level at quasi-linear 
regime (from 4 — 12 /i _1 Mpc). 

While this study focuses on the pairwise velocity of PDF of dark matter, we are in the process to extend the current 
work to study the effect of primordial non-Gaussianity on the pairwise velocity of biased tracers. We are also extending our 
study to examine the effect of primordial non-Gaussianity on the redshift space distortion. The TV-body measurements of the 
pairwise velocity PDF show clearly both the parallel and the perpendicular to the line of separation directions have primordial 
non-Gaussia nity imprints. The effect can be seen in separation as large as r = 50 h~ 1 Mpc. Our analytical model improves 
earlier work I Lam et al. 2010l : Schmidt 2010h - it includes the change in the pairwise velocity PDF in both the parallel and 
the perpendicular directions in the analysis as well as the evolution of the pairwise velocity PDF. As a result it provides a 
much better match to the TV-body measurements compared to the linear theory. This enables our analytical model to improve 
the estimation of the signature of primordial non-Gaussianity in the redshift space distortion. 
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Figure 6. The 1-point PDF of peculiar velocity: the upper panel shows the PDF profile of the peculiar velocity when f n i = 0. The 
open circles show the measurement from TV-body simulation and the solid curve shows the Maxwellian distribution with <r equal to 
the velocity dispersion measured from the simulation. The lower panel shows the ratio of the PDF for / n ; = ±100 to f n i = 0: upward 
pointing triangles are the ratio for /„; = 100; downward pointing for /„; = —100. Particles with large velocities (> lOOOkm/s) are mostly 
populated in massive halos. 
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APPENDIX A: DERIVATION OF MASS WEIGHTED LINEAR PAIRWISE VELOCITY PDF 



This appendix derives the expression of the linear pairwise velocity PDF for /„; 7^ (equation (|20[) '). The following derivation 
also applies to the derivations of the linear velocity difference PDF when /„; 7^ (equation (|13|) . by setting S — S' = in the 
following derivation) and the linear pairwise velocity PDF when f n i = (equation (|17|) ). 

The first s t ep is to generalize the expression for the line-of-sight pairwise velocity generating function M (equation (8) 



Scoccimarrd (|200J)) to the 3d generating function of pairwise velocity: 

Z(X; r) = [1 + ar)]M(\; r) = (exp(A ■ v)[l + 8(x)][l + 8(x')]), 



(Al) 



where A • v = X\\ V\\ + X±„v±„ + X±, v± b and x = x' + r. The characteristics function can be computed using the cumulants 
(see discussion in IScoccimarrd ((200J)): 

Z(i\; r) = exp(exp(iA ■ v)) c [l + (exp(iA ■ v)(S + 5')} c + (exp(iA • v)5) c {exp(iX ■ v)S) c + (exp(iA • v)55') c ] (A2) 



exp 



<«N> 



S3 , < v uO 



3! 



( v \\ v l b ) 



x <l + £(r) + (» u M')(iA||) + 5^ 



'M6 + S'))(i\ j ) + 



(vK5 + 5')) 



(zA,|)(iAiJ 2 + . 



(iX,) 2 + . 



+ 



$Z<«i*)(*Aj) + 



-{ix j y + . 



+ . 



(A3) 



where 8 = 5(a;), 5' = S(x'), and the summation over j denotes summing over all three directions {||, J_ a , .Lb}- The first term 
in the second equity is the linear velocity difference characteristic function while the second term is associated with the mass 
weighting. The characteristic function for the linear velocity difference for /„; = 0, 



.Fo(iA; r) = exp 



(A4) 



whose Fourier transofrm is the linear velocity difference PDF when /„; = (equation (0). Hence we can factorize Tq from Z 
and define Q 

Z(iX;r) = g(iX;r)To(iX;r) (A5) 

So far the expression is exact. Now we take the approximation that only terms up to the first order in /„; is included and 
notice that {v± a S) = (v± b 5) — 0. Hence 



G(iX;r) = 1 + £(r) + «V> + {v [{ 6') + («||5<5')) (iX l{ ) + (u||<5}(«||5'}(iA||) 2 + 



Ml 

3! 



(l + g(r))+ - + 



+ 



<«f> 



<^f> 



.,, (^ || 5> + ( l ; ||( 5'>)(zA | |) 4 + ^(,; ||( 5>( t ; | |<5'>(zA | |) 5 + £ r(zA x J 
where 

r{iX±a) = + (V>L*> + <3|J (1 + CW) 



(»A||) 3 
(A6) 



(iA||)(iA ia ) 2 + M^Ll « V) + (iA||) 2 (iA^) 2 



(A7) 



To get the mass weighted linear PDF we need to do a Fourier transform on Z. While the expression seems very complicated, 
the Fourier transform can be done analytically due to the fact that the Fourier transform of Q is actually a derivative operator 
acting on the uniform weighted linear PDF. Recall that the uniform weighted linear PDF is a product of three univariate 
PDFs, hence the derivative operation can be done easily and the resulting mass weighted linear pairwise velocity PDF is 
equation ([2SJ- 



